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1.  INTRODUCTION 


A  review  of  existing  bottom  backscattering  models  that  are  based  on 
sediment  volume  scattering  was  made  at  the  beginning  of  this  project.  Nolle  et 
al.^  measured  bottom  backscattering  strength  of  water  saturated  sand  in  a  well- 
controlled  laboratory  experiment  and  put  forward  a  bottom  backscattering  model, 
in  which  the  sediment  is  treated  as  an  elastic  and  lossy  liquid.  Using  the 
appropriate  values  for  sound  speed  and  attenuation  in  sand,  it  is  apparent  that 
the  model  predicts  a  sharp  reduction  in  backscattering  strength  below  the  critical 
grazing  angle.  This  feature  is  conspicuously  absent  from  the  corresponding 
laboratory  measurements.  Both  the  theory  and  experimental  results  are  shown 
in  Fig.  1.1.  This  approach  is  essentially  the  same  as  that  of  Jackson^  ®  as  far  as 
scattering  from  smooth  sediments  is  concerned.  More  recently,  Hines^ 
attempted  to  introduce  an  additional  mechanism  that  would  alleviate  the 
discrepancy,  i.e.,  the  coupling  of  acoustic  energy  into  the  sediment  via  the 
evanescent  wave,  also  shown  in  Fig.  1.1.  Doubts  concerning  the  validity  of 
Hines's  formulation  aside,  the  model  still  falls  short  of  agreement  with  the 
laboratory  data. 

It  is  our  hypothesis  that  there  are  two  flaws  in  the  above  models:''"^  (1)  the 
approximation  of  the  sediment  as  a  liquid,  which  breaks  down  at  grazing  angles 
below  the  critical  value,  and  (2)  the  single  scatter  approximation,  which  is  invalid 
where  there  are  numerous  scatterers  in  close  proximity  to  each  other,  in  this  case 
the  sand  grains  themselves.  By  using  the  poroelastic  wave  theory  of  Biot,  and 
the  diffusion  equation,  it  is  our  intention  to  formulate  a  simple,  but  approximate, 
theory  that  can  properly  explain  the  experimental  data. 

The  analysis  for  the  acoustic  backscatter  from  a  smooth  sand  surface 
utilizing  the  diffusion  equation  is  illustrated  in  Fig.  1.2.  An  acoustic  pulse  of 
energy  originating  in  the  water  transmits  across  the  sand/water  interface.  Once 
the  acoustic  pulse  is  within  the  sandy  sediment,  the  energy  diffuses  away  from 
the  pulse,  some  of  which  reaches  the  interface  and  radiates  back  into  the  water. 
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o  EXPERIMENT  (Nolle  et  al.) 

-  THEORY  (Nolle  et  al.,  Jackson  et  al.) 

—  THEORY  (Hines) 
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Figure  1.1 

Review  of  existing  models. 
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wave 


Figure  1.2 

Illustration  of  diffusion  model  of  scattering  in  sandy  sediments. 
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The  assumption  which  therefore  must  be  made  is  that  the  sand  grains  are  in 
contact  with  one  another  to  allow  the  energy  to  be  diffused  through  them.  Once 
the  equation  for  the  diffused  energy  within  the  sand  has  been  established,  the 
energy  must  be  computed  as  it  travels  through  time  and  space  and  the 
necessary  constants  must  be  determined. 
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2.  THE  DIFFUSION  EQUATION 


The  diffusion  equation  has  been  used  in  modeling  ultrasonic  scattering  in 
a  number  of  problems,  including  suspensions,®  polycrystals,®  and  plates/  It  is  an 
approximation  that  is  most  appropriate  for  multiple  scattering  problems  where 
interference  between  scatterers  renders  the  single  scatter  approximation 
untenable.  Due  to  the  high  density  of  scatterers  and  their  close  proximity,  the 
scattering  of  sound  by  sand  grains  in  sandy  ocean  sediments  is  also  expected  to 
be  a  multiple  scattering  problem  that  should  be  modeled  as  a  diffusion  process. 
The  theoretical  derivation  is  as  follows. 

Given  a  density  function  u  of  diffuse  energy,  the  net  energy  flux  E  is  in  the 
opposite  direction  to  the  gradient  of  u, 

E  =  -(Vm)Q?  ,  (2.1) 

where  is  the  diffusion  constant.  In  the  lossless  case,  the  rate  of  change  of  u  is 
then  given  by 


=a'^[v'^u]  .  (2.2) 

at 

In  the  case  of  lossy  diffusion,  there  will  be  a  decay  proportional  to  u  due  to 
absorption,  represented  by  a  decay  coefficient  y^.  In  addition,  there  may  also  be 
an  energy  source  density  term  F.  The  lossy  diffusion  equation  with  a  source 
density  term  is  given  by 


^u  =  cP'(V'^u)-'^u+F  .  (2.3) 

at 

Let  us  consider  two  cases:  the  steady  state  plane  wave  case  and  the  normally 
incident  plane  wave  pulse.  These  two  cases  are  particularly  revealing  of  the 
physical  processes.  The  first  will  show  the  dependence  of  the  scattering  strength 
as  a  function  of  grazing  angle  on  physical  parameters  of  the  medium.  The 
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s©cond  shows  th©  t©mporal  b©havior  of  th©  diffusion  procoss.  Finally,  th© 
practical  matter  of  estimating  th©  vaiues  of  th©  input  parameters  is  considered. 

2.1  STEADY  STATE  PLANE  WAVE  CASE 

In  the  steady  state  case,  there  is  no  change  as  a  function  of  time; 
therefor©  Eq.  (2.3)  reduces  to 

(y\W-y^u+F  =  0  .  (2.4) 

For  a  plane  wav©  propagating  in  the  x-y  plane,  in  which  th©  x  axis  is  horizontal, 
and  the  y  axis  vertical,  there  would  be  no  change  in  the  x  direction.  All  variables 
may  be  expressed  as  functions  of  y.  Given  that  solutions  in  th©  form  of 
exponential  functions  are  known  to  exist,  let  us  express  u  in  terms  of  a  Fourier 
transform. 


•  (2.5) 

Substituting  into  the  steady  state  lossy  diffusion  equation,  an  expression  relating 
Fand  A  is  obtained. 


F  =  j  ([a^k^ +  fyAe^^)dk  • 


The  source  density  F  arises  from  energy  that  diffuses  away  from  the  incident 
wave,  at  a  rate  that  is  proportional  to  the  intensity  of  the  incident  wave.  Since  the 
incident  wave  decreases  exponentially  with  depth,  the  source  density  must  also 
decrease  exponentially.  Let  the  sand/water  interface  be  at  y=0,  values  of  y  less 
than  zero  being  in  the  sand.  Let  us  further  assume  that  the  energy  that  diffuses 
to  the  interface  is  completely  reradiated  into  the  water,  which  implies  that  the 
interface  behaves  as  an  energy  sink.  Then,  source  density  may  be  modeled  as 


(y  < 0) 

-e  ^  (j  >  0) 


(2.7) 
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The  source  density  for  values  of  y  less  than  zero  is  genuine.  The  source  density 
for  y  values  greater  than  zero  is  a  construction  that  effectively  makes  the 
interface  into  an  energy  sink;  the  antisymmetry  ensures  that  the  energy  density  u 
is  always  zero  at  the  interface.  The  source  density  may  also  be  expressed  as  a 
Fourier  transform. 


f  =  r  (Be‘'^)dk 

J—oo 


(2.8) 


Substituting  from  Eq.  (2.8)  into  (2.7),  and  solving  the  integral,  the  Fourier 
coefficient  B  is  given  by 


ik 


(2.9) 


Comparing  Eqs.  (2.6)  and  (2.8),  the  following  relationship  between  the  Fourier 
coefficients  of  the  source  density  and  energy  density  is  found. 


(d‘-k^  +  ‘^)A=B 


(2.10) 


Eliminating  B,  an  expression  for  A  in  terms  of  the  independent  variables  is  found 


ik 


7t{ik  +  P){-ik  +  P)(cp-k^  + 


(2.11) 


Of  particular  interest  is  the  rate  of  diffusion  of  energy  from  the  sand 
towards  the  sand/water  interface.  The  energy  flux  is  related  to  energy  density  by 
Eq.  (2.1).  Substituting  for  A  in  Eq.  (2.5)  from  (2.1 1),  then  for  u  in  Eq.  (2.1)  from 
(2.5),  solving  the  Fourier  integral,  and  then  setting  y=0,  the  energy  flux  towards 
the  sand-water  interface  is  found  to  be 


E  = 


a 

Y+a  p 


(2.12) 
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The  energy  flux  into  the  sand  from  the  sand/water  interface  is  given  by 


(2.13) 


Thus,  the  energy  flux  moving  towards  the  interface  per  unit  energy  flux  moving 
towards  the  interior  of  the  sand  is  given  by 

^  P _  .  (2.14) 

I  7  +a  j3 

Scattering  strength  is  defined  as  scattered  power  per  unit  solid  angle  per  unit 
area  per  unit  incident  intensity.  In  terms  of  E  and  /,  the  bistatic  scattering 
strength  S  is  given  by 


E  sin(0) 

5  =  y0(0,.)(2sin[0J)-^ 

a  p  sin[0/)0  [ai)sin[^o]  ^ 

[y  +  a  p)7i: 

where  ^  (0j)  is  the  plane  wave  intensity  transmission  coefficient  from  water  into  . 
sediment  at  a  grazing  angle  0/;  2sin(0^)  is  the  intensity  transmission 
coefficient  from  sand  to  water  from  a  baffled  aperture  that  is  very  much  smaller 
than  a  wavelength,  at  an  exit  angle  0  o  relative  to  the  horizontal;  0  [0  / )  is 
determined  by  elastic  wave  theory  for  an  elastic  sediment,  or  by  Biot's  theory  for 
porous  sediment  such  as  sand.  Thus,  a  very  simple  expression  for  the  bistatic 
scattering  strength  has  been  obtained  in  terms  of  the  diffusivity  a,  decay 
coefficient  of  the  coherent  wave  as  a  function  of  depth  15,  and  the  absorption 
coefficient  of  the  diffuse  energy  /. 
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2.2.  NORMALLY  INCIDENT  PLANE  WAVE  PULSE 


Consider  a  Gaussian  shaded  pulse  of  width  a  traveling  in  the  negative  y 
direction  at  a  speed  c.  With  reference  to  Eq.  (2.3),  the  source  function  may  be 
expressed  as 


F  =e 


+c  s 


a 


-Pcs 


S[t-s) 


(2.16) 


where  5  is  the  delta  function.  This  expression  was  chosen  to  allow  the  diffusion 
to  be  solved  in  a  stepwise  fashion.  The  first  step  is  to  solve  for  the  solution  of 
each  elemental  time  slice.  Then,  the  elemental  solutions  are  integrated  for  the 
complete  result.  This  method  was  chosen  because  a  simple  analytical  solution  is 
obtainable  for  the  time  slices,  as  follows. 

For  any  particular  value  of  s,  the  above  function  is  nonzero  for  only  an 
instant  when  the  independent  variable  t  is  equal  to  s.  The  contribution  from  that 
instant  us  may  be  treated  as  a  boundary  value  problem  in  which  the  initial  value 
of  Us  is  given.  Since  the  source  term  is  zero  thereafter,  we  may  neglect  the 
source  term  in  the  diffusion  equation,  giving 

[V^u  s]a^-r^u  s  =^u  s  .  (2.17) 


which  has  a  solution  of  the  form 


u  s 


(2.18) 


where  the  separable  component  solutions  are  of  the  form 


r 

r  =<  ®  it  -s  <0)  ^2.20) 

j-it-s)w^  ^0) 

Substituting  from  Eqs.  (2.19)  and  (2.20)  into  (2.18),  and  then  into  (2.17),  the 
following  relationship  is  found  for  the  variables  k  and  w. 

=  .  (2.21) 

The  boundary  condition  is  that,  at  the  instant  when  t-s=0+ ,  the  value  of  Us  is 
simply  equal  to  the  source  function  at  t=s.  Recognizing  that  Eq.  (2.8)  is  a  Fourier 
integral,  the  coefficient  B  may  be  expressed  as  the  inverse  Fourier  transform  of 
the  source  function 


The  solution  is 

-ia'^k^  +  c  (/  k  -P)s 

B  =l££_J _  ■  (2-23) 

Substituting  into  Eq.  (2.19)  and  then  (2.18)  and  solving,  the  result  is 

U  s  = 

(2.24) 

The  total  solution  for  energy  density  is  found  from  integrating  over  all  time  slices 
up  to  the  time  t, 


+r  j  +  a 


f4cr^(— 5  )+( 
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u 


(2.25) 


u  s  d 


s 


Of  particular  interest  is  the  energy  flux  towards  the  water/sediment  interface,  at 
y=0,  which  is  given  by  Eq.  (2.1).  Furthermore,  if  the  interface  is  modeled  as  an 
energy  sink,  we  may  construct  an  antisymmetric  source  in  the  upper  half-space, 
which  has  the  effect  of  doubling  the  energy  flux.  Thus,  substituting  from 
Eq.  (2.25)  into  (2.1)  and  multiplying  by  2,  an  expression  for  the  energy  flux  as  a 
function  of  time  is  obtained. 


E 


at  y=0 


(2.26) 


Substituting  from  Eq.  (2.24)  for  Us,  the  final  solution  is  given  by  the  integral 


E(t)  =  A 


+t]- 


c^s  ^ 


a  c  as  e 


4(x'^[-s  +t]  +  a^ 


c  p  s 


■  d  s 


[4a  +t  \+a 


(2.27) 


This  integral  does  not  have  a  simple  analytical  solution,  and  therefore  it  was 
computed  numerically.  The  integrand  is  monotonic  and  simpie  to  integrate 
accurately. 


If  we  make  the  assumption  that  the  diffusion  constant  is  small  such  that 
A  a^[-s  may  be  approximated  by  a  ^ ,  then  a  simpler  approximate 

solution  is  possible.  Let  Ea  denote  the  approximate  solution 


E,[l)  =  4 


If 


-y2[-S 


case 


■  d  s 


a 


,  (2.28) 
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which  can  be  solved  in  terms  of  error  functions. 


The  total  energy  input  per  unit  area  is  given  by 


noo 

noo 

'  /.oo 

f 

f 

1  F  dt 

dy 

J  ^  oo 

J  —  oo 

> 

(2.29) 


The  integration  with  respect  to  y  should  strictly  be  between  -oo  and  zero.  For 
simplicity,  the  limits  were  extended  to  +  and  -  oo,  so  that  an  error  function  may  be 
avoided.  The  result  is  an  overestimation  of  the  input  energy  by  a  fraction  of 
approximately  afi.  Assuming  that  the  pulse  width  a  is  small  compared  to  the 
decay  constant  1/j3,  the  error  will  be  small.  Substituting  into  Eq.  (2.15),  the 
scattered  power/unit  solid  angle/unit  area/unit  energy  input  is  given  by 


In: 


(2.30) 


2.3  INPUT  PARAMETERS 


The  diffusion  coefficient  is  the  central  parameter  in  this  model  although  it 

does  not  directly  determine  the  scattering  strength.  In  Eq.  (2.15)  it  is  seen  that  in 

2 

the  limit  that  the  diffusion  coefficient  cc  tends  to  infinity,  the  scattering  strength 

saturates  at  an  asymptotic  value,  whereby  all  the  incident  energy  is  scattered 

2 

back  into  the  water.  For  gases,  the  value  of  oc  is  given  by 

a^=c/c/3  ,  (2.31) 

where  d  is  the  mean  free  path,  and  c  the  average  speed  of  particles.  If  we 
assume  that  diffuse  sound  energy  is  transmitted  from  a  sediment  grain  to  its 
neighbors  in  a  similar  manner,  then  the  mean  free  path  d  should  be  equal  to  the 
mean  grain  diameter,  and  c  to  the  speed  of  sound.  Biot's  theory  shows  that 
there  can  be  more  than  one  sound  speed  in  a  porous  sediment.  For  simplicity, 
let  us  make  the  assumption  that  speed  of  energy  transport  can  be  approximated 
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by  a  single  value,  and  that  this  value  is  approximately  equal  to  the  fast  wave 
speed,  which  in  sand  is  typically  in  the  region  of  1700  m/s. 

The  loss  parameter  represents  the  rate  at  which  the  diffuse  acoustic 
energy  is  converted  into  other  forms  of  energy.  This  rate  may  have  some 
relationship  to  the  rate  at  which  energy  is  attenuated  from  the  coherent  acoustic 

wave.  Therefore,  it  is  more  insightful  to  define  7  in  terms  of  a  ratio  ^  y  such 
that 


Ry=y^/{cPn)  ,  (2.32) 

where  c  is  the  speed  in  meters  per  second,  and  P  n  is  the  intensity  attenuation 
coefficient,  in  Nepers  per  meter,  of  a  coherent  plane  wave  traveling  through  the 
sediment:  4.343  P  n  is  equal  to  the  plane  wave  attenuation  in  decibels  per 

meter.  The  value  oi  R  y  must  be  greater  than  or  equal  to  zero.  At  zero,  energy 
is  losslessly  diffused:  at  infinity,  all  diffuse  energy  is  completely  absorbed  by  the 
medium. 

The  attenuation  p  of  the  source  function  is  simply  the  decay  rate  of  the 

plane  wave  intensity  in  the  vertical  direction.  At  normal  incidence  it  is  equal  to 
/5  n  ,  the  attenuation  coefficient  of  the  plane  wave  intensity  in  Nepers  per  meter. 
At  other  angles,  it  is  equal  Xo  P  n  divided  by  the  sine  of  the  transmission  angle, 
measured  relative  to  the  horizontal. 

The  transmission  coefficient  for  sandy  sediments  is  obtained  from  Biot's 
theory.  It  will  not  be  derived  here  because  it  has  been  adequately  covered  in 
readily  available  references.® 
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3.  COMPARISON  WITH  EXPERIMENTAL  DATA 


A  comparison  of  the  model  predicted  backscattering  strength  will  be  made 
with  Lambert's  rule,  and  then  with  backscattering  strengths  measured  by  Nolle 
et  al.  Finally,  the  temporal  nature  of  the  modeled  backscatter  will  be  compared 
with  new  experimental  data  at  normal  incidence. 

3.1  COMPARISON  WITH  EXPERIMENTAL  DATA  FROM  NOLLE  et  al. 

The  experiments  of  Nolle  et  al.  were  done  with  well-graded  sand  samples 
and  narrowband  signals.  Therefore  they  are  directly  comparable  with  the 
theoretical  results  for  discrete  frequencies  and  grain  sizes.  The  attenuation  and 
transmission  coefficients,  P  n  and  0  (0  /),  of  the  coherent  plane  wave  entering 
the  sediment  were  computed  using  the  procedure  described  by  Chotiros^  and  the 
values  in  Table  3.1.  The  parameter  values  used  to  compute  the  diffusion 
constant  a  ^  are  shown  in  Table  3.2.  Little  is  known  about  the  absorption  ratio 
R  Y .  The  diffusion  model  does  not  have  any  means  for  estimating  it  from 
measurable  physical  properties,  aside  from  bounding  its  value  between  zero  and 
infinity.  By  adjusting  the  value  oi  R  y  ,  the  model,  as  expressed  in  Eq.  (2.15), 
could  be  adjusted  to  fit  the  backscattering  strength  measurements.  The  model 
and  experimental  results  at  500  kHz  on  sand  samples  of  four  different  grain  sizes 
are  shown  in  Fig.  3.1.  For  wavelength-to-grain  size  ratios  less  than  10,  the  model 
reaches  saturation,  and  consistently  underestimated  the  measured 
backscattering  strength.  This  does  not  necessarily  mean  that  the  model  is  invalid: 
it  is  likely  that  a  stronger  scattering  mechanism  has  taken  over.  Geometric 
scattering  by  grain  facets,  particularly  those  at  the  water/sediment  interface,  is 
the  most  likely  culprit.  A  similar  set  of  results  for  a  frequency  of  1  MHz  are  shown 
in  Fig.  3.2. 

From  these  results,  an  empirical  relationship  between  the  dimensionless 
absorption  ratio  R  y  and  the  wavelength-to-grain  size  ratio  was  obtained,  as 
shown  in  Fig.  3.3;  in  computing  wavelength,  the  fast  wave  speed  (1700  m/s)  was 
assumed.  Both  the  500  kHz  and  1  MHz  data  appear  to  fall  on  one  continuous 
curve.  The  value  of  /?  y  is  close  to  zero  for  wavelength-to-grain  size 
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Table  3.1 

Water  and  water-saturated  sand  poroelastic  parameters. 


Parameter 

Value 

Units 

Solid  phase  parameters 

kg/m3 

Density 

2650 

Bulk  modulus 

7.00E+09 

Pa 

Fluid  phase  parameters 

kg/m^ 

Density 

1000 

Bulk  modulus 

2.25E+09 

Pa 

Viscosity 

1.00E-03 

kg/m-s 

Frame  parameters 

Porosity 

0.36 

- 

Shear  modulus 

2.61  E+07 

Pa 

Shear  log  decrement 

0.15 

- 

Bulk  modulus 

5.30E+09 

Pa 

Bulk  log  decrement 

0.15 

- 
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Table  3.2 

Diffusion  parameters. 


Parameter 

Symbol 

Value 

Units 

Transport  speed 

c 

1700 

m/s 

Grain  size:  sampie  1 

d 

0.000115 

m 

Grain  size:  sample  2 

d 

0.000164 

m 

Grain  size:  sample  3 

d 

0.00038 

m 

Grain  size:  sample  4 

d 

0.00065 

m 
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Figure  3.1 

Fitting  of  model  backscattering  strength  curves  to  experimental  data 
(Nolle  et  al.)  at  500  kHz  by  adjusting  Ry. 
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Figure  3.2 

Fitting  of  model  backscattering  strength  curves  to  experimental  data 
(Nolle  et  al.)  at  1  MHz  by  adjusting  Ry. 
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ratios  loss  than  10;  at  highar  valuas  of  wavalangth-to-grain  siza  ratio,  ^  y 
appaars  to  incraasa  linaarly. 

3.2  COMPARISON  WITH  LAMBERTS  RULE 

Using  Eq.  (2.15),  tha  paramatar  valuas  in  Tablas  3.1  and  3.2,  and  tha 
ralationship  found  batwaan  R  y  and  wavalangth-to-grain  siza  ratio  shown  in 
Fig.  3.3,  backscattaring  strangth  varsus  grazing  angla  curvas  wara  computad  for 
a  wida  ranga  of  typical  casas.  For  axampla,  tha  scattaring  strangth  curvas  at  1, 
10, 100,  and  1000  kHz  wara  computad  for  a  sandy  sadimant  with  a  grain  siza  of 
0.1  mm.  In  aach  casa,  a  Lambart's  rula  curva  was  comparad  to  tha  diffusion 
modal.  Tha  rasults  ara  shown  in  Fig.  3.4.  In  all  casas,  it  is  apparant  that 
Lambart's  rula  is  a  good  fit  at  grazing  anglas  lass  than  about  30°.  It  is  also  claar 
that  at  10  kHz,  which  corresponds  to  a  wavalangth-to-grain  siza  ratio  of  14782, 
the  backscattaring  strength  is  very  small.  Therefore,  it  may  be  concluded  that  the 
diffusion  model  is  likely  to  be  significant  for  wavelength-to-grain  size  ratios 
between  10  and  10,000.  Above  and  below  these  limits,  it  is  expected  to  be 
overshadowed  by  other  scattering  mechanisms. 

3.3  COMPARISON  OF  DECAY  TIME  WITH  NEW  EXPERIMENTAL 

MEASUREMENTS 

A  unique  feature  of  the  diffusion  model  is  that  it  predicts  a  certain  lag 
between  the  reflected  and  scattered  signals,  which  is  a  function  of  the 
propagation  speed  and  mean-free-path.  While  many  models  can  be  fitted  to  the 
backscattaring  strength  versus  grazing  angle  curves,  few  are  able  to  give  any 
insight  into  the  time  dependence  of  the  scattered  signal.  The  time  dependence  is 
most  clearly  seen  In  the  scattering  of  a  short  plane  wave  pulse  at  normal 
incidence.  An  experiment  was  performed  with  a  wideband  transducer.  A 
narrowbeam  signal  was  used  to  approximate  a  bounded  plane  wave.  Pulse 
compression  filtering  was  applied  to  the  received  signal  to  achieve  an 
approximately  Gaussian  pulse,  with  a  time  constant  of  3.3  ms,  at  a  center 
frequency  of  approximately  150  kHz.  Ten  pings  were  collected  from  different 
parts  of  a  smoothed  sand  surface.  The  sand  had  been  degassed  and  compacted 
by  vibration.  The  coherent  part  of  the  echo  is  the  average  of  the  ten  pings.  The 
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Figure  3.3 

Relationship  between  wavelength-to-grain  size  ratio  and 
the  absorption  ratio  Ry  deduced  from  matching  diffusion  model  of 
backscattering  to  the  experimental  data  of  Nolle  et  al. 
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Figure  3.4 

Example  of  comparison  between  Lambert's 
rule  and  diffusion  theory. 
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scattered  component  is  computed  from  the  rms  difference  between  the  individual 
pings  and  the  average. 

This  experiment  is  different  from  that  of  Nolle  et  al.  in  two  ways:  (1)  a 
broadband  signal  was  used  in  order  to  achieve  a  short  signal  pulse,  and  (2)  the 
sand  was  not  graded  and  contained  grains  of  various  sizes.  Assuming  that  the 
grain  and  wavelength  effects  add  in  a  mean-square  sense,  the  response  of  the 
sand  to  a  broadband  signal  may  be  estimated  by  convolving  E  )  with  the 
density  functions  of  the  grain  size  P  q[(l]  and  signal  spectrum  P  f{f  ],  where 
<?  is  grain  size  in  (j),  and/  is  frequency  in  Hertz.  The  final  result  c  )  is 
then  given  by  the  following  convolution  integral. 


[q]E  [t]d  q 


Pf{f)df 


(3.1) 


The  cumulative  distribution  function  of  the  grain  size  was  measured  at  intervals  of 
1  (j)  and  a  normal  distribution  function  was  fitted  to  it,  as  shown  in  Fig.  3.5,  from 
which  p  q[q]  was  calculated.  The  acoustic  signal  was  digitally  filtered  and  the 
resulting  spectral  density  function  of  signal  spectrum  p  f{f  ]  as  shown  in 
Fig.  3.6.  These  functions  were  used  to  numerically  compute  Eq.  (3.1),  in 
conjunction  with  the  parameter  values  in  Tables  3.1  and  3.2,  and  values  oi  Ry 
as  a  function  of  wavelengths-to-grain  size  ratio  as  shown  in  Fig.  3.4. 

The  theoretical  and  experimental  results  are  shown  in  Fig.  3.7.  The  signal 
envelopes  are  plotted  on  a  decibel  scale.  The  theoretical  and  experimental 
results  appear  to  be  consistent  for  the  most  part.  The  disagreement  is  mostly  in 
the  vicinity  of  t=0,  and  this  is  likely  caused  by  an  edge  effect  that  has  been 
ignored  in  the  theoretical  model;  the  model  assumes  that  the  pulse  is  a  fully 
formed  Gaussian  shaded  function,  as  stated  in  Eq.  (2.16).  In  reality,  at  t=0  only 
half  the  pulse  has  entered  the  sediment,  and  there  are  likely  to  be  further 
complications  due  to  surface  effects.  The  model  appears  to  overestimate  the 
scattered  signal  in  the  vicinity  of  /=0.  The  apparent  peak  in  the  experimental 
measurements  of  the  random  component  of  backscatter  about  t=0  should  be 
ignored  because  it  is  most  likely  caused  by  surface  scattering,  due  to 
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Figure  3.5 

Measured  cumulative  distribution  of  grain  sizes  and  best  fit 
normal  probability  distribution  (mean  =  ^.7^,  standard  deviation  =  0.75(1)). 
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Relative  signal  level  -  dB 


Figure  3.7 

Comparison  of  diffusion  model  with  random  part  of  measured 
reflected  signal  at  normal  incidence  from  a  smooth  sand  surface. 
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residual  surface  roughness.  The  experimental  data  are  only  valid  up  to 
f  =  100  \is,  because  multipaths  within  the  water  column  begin  to  appear  at  this 
point.  Generally,  the  model  and  experimental  data  are  in  agreement. 
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4.  CONCLUSIONS 


A  scattering  model  from  water-saturated  sand  with  a  smooth  surface, 
based  on  the  diffusion  theory,  was  developed.  It  is  applicable  to  acoustic 
scattering  from  sandy  ocean  sediments.  It  has  the  following  features. 

(1)  Scattering  strength  versus  grazing  angle  curve  closely  resembles 
Lambert's  rule. 

(2)  Its  predictions  can  be  closely  matched  to  laboratory  measurements  by 
Nolle  et  al. 

(3)  The  comparison  with  experimental  data  yielded  a  tentative  relationship 
between  the  absorption  and  the  wavelength-to-grain  size  ratios,  which 
indicated  that  scattering  by  diffusion  is  expected  to  be  significant  for 
wavelength-to-grain  size  ratios  between  10  and  10,000. 

(4)  The  temporal  characteristics  of  the  scattered  signal  are  consistent  with 
new  experimental  data,  in  which  the  random  part  of  the  signal  was 
detected  and  found  to  be  consistent  with  the  model  prediction  in  both 
amplitude  and  decay  constant. 

(5)  Since  a  large  proportion  of  shallow  water  regions  of  the  globe  are  sandy, 
this  model  is  expected  to  be  a  significant  component  in  a  comprehensive 
bottom  scattering  model  for  modeling  high  frequency  sonar  systems. 

This  research  project  is  closely  related  to  another  basic  research  project 
sponsored  by  ONR  Code  321 OA,  to  extend  Biot's  theory  for  the  purpose  of 
modeling  scattering  by  sediment  grains.  It  is  expected  that  when  the  ONR 
sponsored  project  is  completed,  it  will  supply  a  sound  basis  for  computation  of 
the  input  parameters  for  the  diffusion  model. 

The  model  for  multiple  scattering  by  sediment  grains  developed  in  this 
report  will  be  transitioned  into  the  larger  bottom  scattering  model  currently  being 
assembled,  under  the  MTEDS  project.  The  total  model  contains  components 
developed  in  past  basic  research  projects  sponsored  by  the  NRL  High  Frequency 
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Acoustics  Program  over  the  past  five  years,  including  the  model  for  scattering  by 
sediment  gas  bubbles.  The  final  objective  is  to  obtain  OAML  certification. 
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